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摘 要 :针对 非 恭 浅 阻 尼 结 构 基于 Kanai-Tajimi 谱 的 卷 积 -微分 混合 动力 方程 解法 较 繁琐 的 问题 , 提 
出 了 一 种 新 的 简明 封闭 解法 。 非 和 蒜 滞 阻 尼 模 型 能 较 好 地 模拟 实际 工程 材料 的 阻尼 特征 , 常 以 指数 
型 核 函 数 的 卷 积 形式 表示 ,给 出 其 对 应 的 微分 型 本 构 关系 。Kanai-Tajimi 谱 随 机 地 震动 模型 能 较 好 
地 描述 场地 的 随机 地 震动 特性 ,其 工程 应 用 时 所 获得 的 结构 地 震动 响应 表达 式 复杂 ,但 其 可 用 滤波 
方程 表示 为 基于 白 噪声 激励 的 随机 过 程 。 利 用 非 锻 滞 阻尼 结构 的 微分 型 本 构 关系 和 Kanai-Tajimi 
谱 的 滤波 方程 ,将 基于 复杂 地 震动 激励 卷 积 -微分 型 动力 方程 转化 为 基于 白 噪 声 激励 的 全 微分 型 动 
力 方程 组 ;然后 基于 复 模 态 法 及 Dirac 函数 的 性 质 ,获得 了 耗 能 结构 系列 响应 (位 移 及 速度 )0 ~2 阶 
谱 矩 的 简明 解 ; 最 后 ,分 析 了 基于 首 超 破坏 准则 及 Markov 假定 的 非 黏 滞 阻 尼 结 构 的 动力 可 靠 度 。 
对 一 算 例 运用 所 提 方 法 与 虚拟 激励 法 进行 对 比分 析 , 证 明 所 提 封 闭 解 的 正确 性 和 高 效 性 ,并 可 作为 
虚拟 激励 法 精度 的 验证 方法 。 
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Abstract: In order to resolve the tedious problem of convolutional-differential hybrid equation solution 
based on Kanai-Tajimi spectrum of non-viscous damping structures ,a new concise closed solution method 
is proposed. The non-viscous damping model usually expressed in the convolution form of an exponential 
kernel function can well simulate the damping characteristics of practical engineering materials ,while the 
corresponding differential constitutive relationship of the damper is given in this paper. The random 
ground motion characteristics of the site can be better described by Kanai-Tajimi spectrum random ground 


motion model ,from which the structural ground motion response expression obtained in its engineering ap- 
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plication is complex ,but it can be expressed as a random process of white noise excitation by filtering e- 


quation. Using the differential constitutive relationship of the non-viscous damping structures and the filte- 


ring equation of Kanai-Tajimi spectrum ,the complex ground motion excitation convolution-differential dy- 


namic equation is converted into a fully differential dynamic equation system excited by white noise , and 


then the complex modal method and the concise characteristics of the white noise excitation are used to ob- 


tain the variance of the structural displacement and the concise closed solution of 0-2 order spectral mo- 


ment. Finally ,the dynamic reliability of the structure is analyzed based on the first exceedance failure cri- 


terion and Markov distribution assumptions. A comparative analysis of the method in this paper and the 


Virtual excitation method for a calculation example shows that the closed solution proposed in this paper is 


correct and efficient ,and can be used as a verification method for the accuracy of the virtual excitation 


method . 


Key words: Kanai-Tajimi spectrum;non-viscous damping structure ;filtering equation ; concise closed solu- 
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忌 基于 地 震 灾害 所 导致 的 人 员 伤亡 和 财产 损失 的 
开 重 性 ,防震 减灾 "研究 成 为 当下 重要 的 研究 课 
题 () 阻 尼 是 振动 控制 领域 中 用 来 描述 振动 过 程 能 量 
耗 艇 的 关键 参数 , 对 结构 动力 响应 的 影响 至 关 重 
要 学 。 但 鉴于 阻尼 的 作用 机 理 异 常 复杂 , 目前 仍 没 
有 统一 的 阻尼 模型 。 工 程 上 应 用 较为 广泛 的 是 茜 滞 
陋 忆 模型 ,该 模型 认为 阻尼 力 与 速度 成 正比 ,未 考虑 
握 吉 滞后 现象 , 与 实际 不 符 。Biot 根据 阻尼 力 与 
质 诅 速度 的 时 间 历 程 相关 的 特性 , 提出 了 一 种 用 质 
点 氢 动 速度 与 核 函 数 ( 常 为 指数 型 ) 卷 积 的 阻尼 力 
模 表 。Liu!” 研究 了 时 程 地 震动 激励 下 非 符 滞 阻尼 
结 物 的 动力 响应 分 析 。 段 忠 东 等 中 利用 拉 普 拉 斯 变 


域 法 和 时 域 法 “” 。 在 频 域 法 中 , 虚拟 激励 法 “” 
和 传递 函数 法 “为 其 主要 方法 ,但 上 述 方法 分 析 
结构 响应 方差 或 0 ~2 阶 谱 矩 时 需要 对 结构 响应 的 
功率 谱 在 频 域 [10, + %) 进行 数值 积分 , 虽 提出 了 一 
些 快 速算 法 ” ,但 对 于 复杂 的 激励 模型 或 超大 型 
结构 的 动力 响应 分 析 仍 然 存 在 计算 效率 低 及 精度 不 
足 的 问题 。 在 时 域 法 中 , 复 模 态 法 是 典型 方 
法 ”| ,其 将 任意 形式 的 线性 振动 体系 进行 复 模 
态 解 耦 ,将 结构 的 随机 响应 的 协 方差 表示 成 脉冲 函 
数 与 激励 协 方差 的 二 重 积分 , 故 其 应 用 的 前 提 是 激 
励 要 有 协 方差 且 越 简单 越 容 易 求 封 闭 解 。 葛 新 广 
等 5220 分 别 利用 Kanai-Tajimi 谱 和 Clough-Penzien 


换 研 究 了 非 纤 滞 阻 尼 模 型 中 的 阻尼 系数 识别 问题 ， 
研究 表明 非 条 灌 阻尼 模型 能 较 好 地 反映 实际 结构 的 
耗 角 性 能 。 文 献 [9] 针 对 非 纤 滞 阻 尼 结 构 的 地 震动 
时 程 激励 ,提出 了 一 种 指数 核 函 数 和 高 斯 核 函数 结 
构 系统 直接 积分 的 简化 方法 。 以 上 研究 都 是 针对 非 
符 灌 阻尼 结构 基于 时 程 激励 下 的 动力 响应 方法 进行 
的 研究 。 
通过 对 大 量 的 地 震动 观测 资料 研究 表明 ,地 震 
动 具有 显著 的 随机 性 ,并 且 各 种 研究 还 提出 了 数 种 
地 震动 随机 模型 "2 。Kanai-Tajimi 模型 “001 能 考 
虑 场地 的 动力 特征 对 地 震动 的 影响 ,有 着 广泛 的 工 
程 应 用 。 目 前 基于 Kanai-Tajimi 谱 激 励 下 的 结构 地 
震动 响应 表达 式 均 较为 复杂 ,不 利于 工程 应 用 。 故 
研究 基于 Kanai-Tajimi 模型 激励 下 非 蒜 滞 阻 尼 结 构 
简明 响应 解 ,具有 重要 的 工程 实践 意义 。 

研究 随机 激励 引起 结构 地 震动 响应 的 方法 有 频 


谱 的 滤波 方程 ,将 复杂 的 地 震动 激励 转化 为 具有 
Dirac 函 数 表 示 的 协 方差 , 给 出 了 时 域 法 求解 随机 响 
应 的 新 思路 。 

针对 以 上 非 笑 滞 阻 尼 结 构 地 震动 响应 研究 的 不 
足 , 本 文 利用 文献 [19-20 ] 所 提 方法 研究 了 指数 型 核 
函数 的 非 医 滞 阻 尼 结 构 基 于 Kanai-Tajimi 谱 随机 地 
震动 响应 的 简明 封闭 解 。 首 先 利 用 Kanai-Tajimi 谱 
的 滤波 振动 方程 将 复杂 的 地 振动 精确 地 表示 成 基于 
Dirac 函数 表示 的 白 噪 声 的 随机 地 震动 ;其 次 将 指数 
型 非 黏 滞 阻 尼 结 构 积 分 型 的 本 构 关 系 转换 成 等 效 的 
微分 型 本 构 关系 ;然后 联 立 地 震动 的 滤波 振动 方程 
及 微分 型 的 非 黏 汪 阻尼 本 构 关 系 ,获得 结构 地 震动 
基于 白 噪 声 激励 的 全 微分 方程 组 ,并 运用 复 模 态 法 
进行 解 耦 ,从 而 获得 结构 位 移 速度 等 响应 的 方差 、 功 
率 谱 和 0 ~ 2 阶 谱 矩 的 简明 解 ; 最 后 研究 基于 首次 超 
越 破坏 准则 及 Markov 超越 分 布 规则 的 动力 可 
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1 基于 Kanai-Tajimi 谱 韭 锋 滞 阻尼 
结构 地 震动 方程 重 构 
对 于 单 自由 度 非 昔 滞 阻 尼 结 构 , 地 震动 方程 为 
mi + [CC DEC dr + he =— mi (1) 
式 中 :mm 分别 为 结构 的 质量 和 刚度 ;x 、% x 分 别 为 
结构 相对 于 地 面 的 加 速度 、 速 度 、 位 移 ;x*, 为 Kanai- 
Tajimi 谱 地 面 运动 加 速度 ;G(?) 为 非 秋 浪 阻 尼 系 统 
的 阻尼 核 函 数 ,一般 采用 指数 形式 , 即 
G(t) = 人 6G.(1) = > coe™ (2) 


式 电 :c 为 阻尼 常数 ;n 为 阻尼 作用 下 松弛 因子 的 数 
量 43 为 松弛 因子 。 


t 


加 令 P = case “4(r)dr ,并 对 已 求 导 得 


9 忆 +awP -cax(bi) =0 (3) 
CODKanai-Tajimi 谱 地 震动 模型 可 将 地 面 运动 表示 
关 基 于 和 白 噪声 激励 的 随机 过 程 , 即 
人 % = 也。 十 LER (4) 
| Ls 2E OU + Os = -in (5) 
式 香 ,is wu 分 别 为 场地 土 相对 于 基 岩 运动 的 加 
速 并 .速度 .位移 ; zx 为 基 岩 的 绝对 加 速度 , 且 具 有 
白 碟 声 特征 ;#, .w, 分 别 为 基 岩 以 上 场地 土 的 阻尼 比 
和 左 率 。 
联 立 式 (2) ~ 式 (5) , 式 (1) 用 矩阵 形式 表示 , 即 


My +Ky=run (6) 
y=[% us % us PP, P ] (7) 
式 中 
r=[-1 000 0 … OJ 44) x13 
[0 10 2éw, 0 0 01 
m 0 0 -2nco 0 0 0 
0 0 1 0 0 0 0 
M=|0 0 0 1 0 0 0 
0 0 0 0 1 0 0 
0 0 0 0 0 
[0 0 0 0 0 和 … 1 Jo x Cn+4) 


人 期 刊 sl 


0 0 0 mw 0 0 0 
0 0 -no 1 0 0 
-1 00 0 .000 
K=| 0 -10 0 00 0 
-Wa “0 0 a 0 
: 00 0 0 
-ca, 0 0 0 0 a 


Ss: 泊 ni (n+4) x (n+4) 


2 复 模 态 解 耦 


运用 复 模 态 理论 ”” , 式 (6) 存 在 特征 值 矩 阵 
P. 右 特征 向 量 矩 了 泗 U 和 左 特征 向 量 和 矩阵 站 ,使 其 解 
耦 , 即 


P=(VMU) 一 (KUD) (8) 
式 中 了 为 对 角 阵 。 
引入 复 模 态 变换 
y=Uz (9) 


式 中 z 为 复 模 态 变量 。 
将 式 (9) 代 入 式 (6) ,并 在 其 前 左 乘 Y 得 


VMUzZ+V'KUz=Vriun (10) 
由 式 (8) 可 知 , 式 (10) 可 改写 为 
z+Pz=n un (11) 


式 中 y= (VMU) ”+ (Vr)。 
由 于 己 为 对 角 和 矩阵 ,根据 哈 梅 积分 , 式 (11) 的 
分 量 为 


2 neinct — 7)d7 (12) 
0 


式 中 Zi、\ Ni 分 别 为 Z,\n 的 分 量 。 
根据 式 (7)、 式 (9) 及 式 (12), 则 结构 相对 于 地 
面 的 位 移 及 速度 可 分 别 表示 为 


n+4 ‘ 


“= Wt = DAssleriia(t -7)dr (13) 
i=1 0 


n+4 t 


w= 本 区 三 ijeia(t -7)dr (14) 
z=] 0 


式 中 :i 为 地 震 激励 作用 时 间 ;u; 为 右 特 征 向 量 矩 阵 
U0 的 第 i 行 向 量 ;响应 的 模 态 强度 系数 和 ,为 

Aij = Uj; * 7; (15) 
式 中 尺 , ,为 右 特征 问 量 矩 阵 上 U 第 i 行 、 第 j 列 的 元 素 。 


3 ”结构 响应 的 方差 及 功率 谱 


观察 式 (13) 和 式 (14) ,结构 位 移 和 速度 统一 表 


示 为 
y(1) = Dp jai -Dar = > 
| - 
(16) 
由 式 (16) 可 得 ,结构 的 响应 分 量 为 
站 (人 人 _7)d7 (17) 


式 中 :=A 为 结构 位 移 的 模 态 强度 系数 ;p= 入 1， 
为 结构 的 速度 的 模 态 强度 系数 ;7 为 时 间 差 。 
结构 响应 的 协 方 差 为 
Cy(7) = EL[Y(1) Yr + 7)] 
= 5S 5 ELY,() YG + (18) 
由 式 (18) 可 得 ,结构 响应 分 量 的 协 方 差 为 
BY (t)Y,(rt + 7)] 


Hon 人 | ee X 五 [xn(t 一 wun (1 +7 -wv) ldudv 
> (19) 
一 恨 据 相关 函数 与 协 方差 的 关系 , 式 (19 ) 改 写 为 
C i (t)Y,(r + 7)] 
> 4 Ki | | eerCicu 十 T 一 v) dudv (20) 
二 为 积分 变量 。 
- 白 品 声 随机 激励 的 协 方差 C ;(7) 可 表示 为 
Ci(T) =2mS05(7) (21) 
8， So 为 地 震动 强度 常数 ; ;8(7) 为 Dirac 函数 。 
殊 将 式 (21) 代 入 式 (20) 得 
CELY,(1)Y,(t + 7)!] 
a | | eer” x 6(u+7T-v)dudv (22) 
根据 Dirac 函数 的 性 质 , 式 (22 ) 化 为 一 重 积 分 
ELY,(D) YC + 7)1 = 255o pup er ed 


(23) 
对 式 (23 ) 积分 得 
ELY, (Y(t +7)]1= 2T9o nn: a (24) 
由 式 (18) 式 (24) , 单 自由 度 非 竺 ; 篆 结 构 基于 
Kanai-Tajimi 随机 激励 下 的 响应 表达 式 为 
Cy(7) = ELY, CE + 7)] 
VT Leb pr 

-2 (5) 

本 SN 人 LE Ni 
人 1 Py TDPi (26) 


式 (25 ) 可 表示 为 
Cy(T7) = -2TSu > Die™ (27) 
当 7=0 时 ,结构 响应 的 协 方差 等 于 结构 啊 应 的 


方差 , 即 
oy” = Cy(0) =-— 215, 5 D, (28) 
式 (28) 可 表示 结构 振动 特征 值 指数 函数 的 线 
性 组 合 
由 en Khinchin 关系 ,结构 响应 的 单 边 功率 
谱 Sy(o) 为 


Sy(w) = | Cr)eos(wr)dr (29) 
将 式 (27) 代 入 式 (29 ) 得 
Sy(w) = 250 > -2 (30 ) 


w” +p,” 
由 式 (16) 式 (28) 、 式 (30) 可知 ,结构 基于 Ka- 
nai-Tajimi 激励 下 的 位 移 及 速度 的 方差 和 功率 谱 的 
表达 式 简 洁 明 了 。 


4 结构 响应 的 0~2 阶 谱 矩 分 析 


由 随机 振动 理论 二 ,结构 响应 Y(t) 的 i 阶 谱 矩 


Bi 可 表示 为 
B; = 2 | wiSyCw) do (31) 
将 式 (30) 代入 式 (31) ,并 令 i=0, 则 结构 响应 
的 0 阶 谱 矩 可 表示 为 
加 > pr 
4 | 2 ee de (32) 
对 式 (32 ) 积分 可 得 
Bo = 25s, 5D, (33) 


平稳 激励 下 的 结构 响应 总 的 0 阶 谱 矩 B。 , 通 
比较 式 (33 ) 及 式 (28 ) , 两 式 完全 相同 ,说明 
we 阶 谱 矩 的 正确 性 。 由 随机 振动 
理论 '" ,结构 响应 的 速度 的 0 阶 谱 矩 等 于 其 位 移 的 
2 阶 谱 矩 , 即 


B, =o7y(0) (34) 
式 中 了 为 了 对 时 间 + 的 一 阶 导数 。 
由 谱 和 矩 的 定义 可 得 ,1 阶 谱 矩 B, 可 表示 为 
加 归 让 Pi 
Bi = 45, > | 7 Bd (35) 
对 式 (35) 进行 积分 得 


n+4 
B! =- 25o > Dipilnps + 


lim 


250 一 二 zine + p?) F Dw, (36) 


由 文献 [12] 可 知 , 实 Dip, =0, 且 实际 频率 不 

可 能 无 限 大 , 故 1 阶 谱 和 矩 8 可 表示 为 
B =- 25, 5 Dpilnp? (37) 
由 起 (33) 起 (34) 式 (37) 可 知 ,本 文 所 获得 的 
结构 基于 Kanai-Tajimi 激励 下 的 相对 于 地 面 位 移 的 


0 ~2 阶 谱 和 矩 均 有 闭 式 解 ,表达 式 简 洁 明 了 ,与 文献 
[12-13 ] 所 提 方 法 相 比 ,本 文 方法 不 需要 积分 。 


5 :动力 可 靠 度 


工程 上 ， 线性 结构 可 靠 度 的 分 析 主 要 是 基于 首 
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s[m,al =125s  ， 


数 G(1) =c(ae “+Qe 2 ) ,其 阻尼 核 函 
值 参考 文献 [6] ,c=9.5x1l0”N . 
oo =255 s 。 
6.1 地 面 加 速度 %, 的 功率 谱 对 比分 析 

地 面 加 速度 可 由 式 (4 \ 式 (5) 式 (7) 式 (16) 
联 立 可 得 , 即 


本 2 
WE 26 ,eus WU = 


n+4 " 
>» je — T)d7 
i=1 0 


(41) 
式 中 «; 为 地 面 加 速度 的 模 态 强度 系数 , 且 k; = 
-25.oA -wzA4 ,其 中 和 Ai 的 含义 见 式 (15) 。 
按 本 文 方法 ,地面 加 速度 模 态 强度 系数 k; 如 表 

1 所 示 , 其 中 等 效 阻尼 比 &, =real(p;)/|p;|。 
由 表 1 可 知 , 模 态 1 的 特征 值 为 复 共 斩 项 ,其 等 
效 频率 为 w, ,等 效 阻尼 为 &,, 且 x, 的 模 态 强度 系数 


超 破坏 准则 及 Markov 分 布 规则 ,其 中 文献 [23] 
提出 了 双 侧 相同 超越 界限 的 可 靠 度 计算 公式 


人 (38) 

E | | ld) | 
Nop -= | > 2 
1 -exp| =| 

(39) 

式 帅 ， R(n) 为 第 nn 年 的 界限 值 ;7 为 一 次 激励 的 持续 


时 闻 ;o, .0 分 别 为 响应 * 和 其 变化 率 的 方差 ;9, 为 


请 大 数 ,其 计算 公式 为 
_/_h 
A BoB, 
式 中 Bo、B1、B, 分 别 由 式 (33) 、 式 (37)、 式 (34) 计算 
得 出 。 
6 算 例 


建立 在 7 度 区 二 类 场地 上 的 单 层 结构 , 梁 \、 柱 尺 
寸 分 别 为 300 mm x550 mm .600 mm x 600 mm , 层 高 
为 6m, 跨 度 宽 为 6m, 按 刚性 楼 板 层 假定 。 结 构 质 量 

m =1.5 x105 kg, 结 构 抗 侧 刚 度 上 =3.6 x10 N/m, 地 
震动 激励 为 Kanai-Tajimi 谱 , 其 参数 取 值 :场地 阻尼 
比 上 =0. 72 ,特征 频率 w. =15. 71 rad/s , 白 噪 声 强 度 
因子 $6 =31.76 x10… m/s ; 非 秋 滞 阻尼 模型 核 函 


不 等 于 零 , 而 其 余 模 态 的 强度 系数 均等 于 0。 这 与 
场地 的 地 震动 仅 与 场地 的 动力 特征 有 关 、 而 与 结构 
的 特征 无 关 所 决定 ,从 侧面 论证 了 本 文 方法 的 正 
确 性 。 


表 1 地 面 加 速度 %, 的 模 态 强度 系数 
Tab.1 Modal intensity coefficients of ground acceleration %, 
等 效 ”等 效 
模 态 特征 值 _ 模 态 强度 
网 频率 ”阻尼 比 和 
|pi | 他 


1 -11.3111 干 10.9023i 15.71 0.72 -11.3112 干 0.416 Si 


2 —0.6169+0.1596i1 15.98 0.04 0 
3 —124.3675 124.36 1.00 0 
4 —254.3691 254.36 1.00 0 


本 文 方法 中 ,基于 Kanai-Tajimi 谱 的 地 面 绝对 加 
速度 激励 功率 谱 ， 0) . 式 (41) 可 得 


6 
KiKi Fk 
Sa ) = 249 (42) 
w eg 


传统 方法 中 ,基于 Kanai-Tajimi 谱 的 地 面 绝对 加 
玉 度 激励 功率 诺 其 表达 式 为 
Ss wt +4é,0.0 (43) 


(oz 一 oO) +48zaw? 2 

图 1 旨 会 制 了 式 (42) 所 表示 的 本 文 方法 和 式 
(43) 所 表示 的 传统 方法 的 绝对 加 速度 功率 谱 密度 函 
数 图 ,绘图 区 间 取 [0,60]rad/s, 绘 图 间距 为 2 rad/s。 
观察 图 1 ,两 者 完全 重合 ,验证 了 本 文 方法 的 分 析 功 
率 谱 的 正确 性 。 
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一 一 相交 方法 
一 一 传统 方法 


rad ') 


Soy(lO mm + s+ 


0 10 20 30 40 50 60 
witrad vs 
色 1 地 面 加 速度 功率 谱 对 比 图 
Fig. 1 Power spectrum comparison of ground acceleration 
6.2 结构 位 移 功 率 谱 对 比分 析 
本 文 方法 中 ,基于 Kanai-Tajimi 谱 的 单 自由 度 结 
构 的 相对 于 地 面 的 位 移 功率 谱 , 由 式 (30) 得 
(© Se WH: Pr 
S,(w) = 2 2 这 a (44) 
司 传 统 方法 中 ,系统 的 幅 频 特性 及 其 基于 Kanai- 
Tiihi 谱 的 单 自 由 度 结构 的 相对 于 地 面 的 位 移 功 率 
谱 贰 表达 式 分 别 为 


1 1 8 和 
9 (w) = Oem (45) 
(0w) H(jo)H( -jw) 
式 中 H (jw) = (km + + ol 入 
SS Ql +]w Q2 + 
ST. 
: 12 
7 一 9 一 本 交 方 法 
可 一 # 一 从 统 方法 


SMIDO mm +s* rad 


0 10 20 30 4d0 50 o0 
mitrad 。s 


图 2 结构 位 移 功率 谱 对 比 图 

Fig.2 Comparison of PDF of structural displacement 

图 2 为 本 文 方法 与 传统 方法 中 结构 位 移 的 功率 
谱 密 度 函 数 图 。 由 图 2 可 知 ,本 文 方法 与 传统 方法 
的 曲线 完全 吻合 ,验证 了 本 文 方法 计算 结构 响应 功 
率 谱 密 度 函 数 的 正确 性 。 本 文 方法 为 式 (44) , 相 较 
传统 方法 的 式 (45 ) 更 为 简明 。 
6.3 计算 效率 及 精度 

根据 虚拟 激励 法 ,结构 响应 的 0 ~ 2 阶 谱 矩 计 算 


式 为 
Bs = oz = 人 S.C) dw (46a) 
Bb! = | S,(w) wdw (46b) 
康吉 品 三 | SCo)azda (46c) 
式 中 S,(w) 见 式 (45)。 


观察 式 (46) 的 响应 值 计 算 , 仅 能 采用 数值 方法 
在 [0, +%) 进行 求解 ,显然 是 无 法 实现 的 。 由 图 2 
可 知 , 随 着 频率 w 的 增 大 ,功率 谱 越 来 越 小 。 本 算 例 
中 ,积分 区 间 取 [0,250] rad/s, 积 分 步 长 取 3 种 情况 
对 比 验证 虚拟 激励 法 的 精度 问题 :中 积分 步 长 为 
1 rad/s; 四 积分 步 长 为 0.5 rad/s; 加 积分 步 长 为 
0.1 rad/s。 表 2 列 出 两 种 方法 的 0 ~2 阶 谱 和 矩 精度 对 
比 , 表 3 列 出 两 种 方法 的 0 ~2 阶 谱 和 矩 计算 效率 对 
比 。 


表 2 0 ~2 阶 谱 算 精度 对 比 


Tab.2 Accuracy comparison of 0-2 order spectral moments 


方法 Bo/m” Bi/(m*s™') Bs/(m s 2) 
本 文 方法 ”4.53016x10-5 7.00194x10-4 0.011 13 
OO 4.70009 x10-5 7.27642x10-4 0.011 57 
@ 4.53303x10-5 7.00661 x10-4 0.011 14 
®@ 4.53016x10-5 7.00194x10-? 0.011 13 


由 于 虚拟 激励 法 的 谱 矩 计算 采用 数值 积分 ,其 
分 析 精 度 随 着 积分 步 长 减 小 而 提高 。 由 表 2 可 知 ， 
虚拟 激励 法 的 0 ~2 阶 谱 矩 计算 精度 , 随 着 积分 步 长 
从 1 rad/s、0.5 rad/s 和 0. 1 rad/s 的 减 小 ,所 得 谱 矩 
值 逐 渐 接 近 本 文 方法 所 得 谱 算 值 。 理 论 上 说 明 本 文 
方法 为 封闭 解 , 且 精 度 较 好 。 
表 3 0 ~2 阶 谱 矩 计算 效率 对 比 


Tab.3 Comparison of computing efficiency 


of 0 ~2 order spectral moments 


方法 虚拟 激励 法 本 文 方法 
OO 0.087 s 
©® 0.101s 0.024s 
® 0.229 s 


由 表 3 可 知 ,虚拟 激 励 法 计算 谱 矩 的 耗 时 随 着 
了 分 步 长 的 减 小 而 增加 。 由 表 2 可 知 , 积 分 步 长 为 
1 rad/s 时 ,精度 较 差 ,积分 时 间 是 本 文 方法 的 3.6 
音 ;积分 步 长 为 0.1 rad/s 时 在 取 6 位 有 效 数 值 的 情 
况 下 ,两 者 完全 一 致 ,而 耗 时 却 是 本 文 方法 的 9.5 


音 , 因 此 本 文 方法 的 计算 效率 较 高 。 
6.4 结构 动力 可 靠 度 

结构 的 层 间 位 移 取 层 高 卫 的 1/250,H=6.0 m， 
按 上 下 界限 均 取 24 mm ,荷载 持续 时 间 了 = 15 s, 可 
徘 度 分 析 的 基本 参数 如 表 4 所 示 。 

由 于 侧 移 方差 为 6.61 m, 远 小 于 24 mm, 但 可 靠 
度 也 只 有 0.947 35 ,不 满足 现行 规范 要 求 , 说 明 动 力 
对 结构 的 影响 不 能 仅 根据 侧 移 值 来 判定 。 

表 4 响应 方差 及 可 靠 度 计算 参数 


Tab.4 Response variance and reliability 


calculation parameters 


Oi/m 


oi/(m:s-!) gs pr 


6.613 52 x10 -3 7.82515 x107™? 0. 167 02 0.947 35 


构 基 于 Kanai-Tajimi 谱 随 机 地 震动 响应 进行 了 研究 ， 
结论 如 下 。 
1) 非 秋 沾 阻尼 模型 采用 指数 型 核 函数 时 ,推导 
无 虚 精 确 的 一 阶 微分 型 本 构 关系 ,便于 与 结构 体系 
的 动力 方程 联合 求解 结构 响应 。 
2) Kanai-Tajimi 谱 的 滤波 方程 将 复杂 的 地 面 运 
动 激励 精确 地 表示 成 具有 简明 协 方差 关系 的 白 噪 声 
随 丢 地 震动 ,其 与 结构 的 运动 方程 联合 建 模 ,利用 复 
模 态 方法 从 时 域 角度 便于 获得 结构 响应 方差 的 简明 
闭 屋 解 和 从 频 域 角度 获得 结构 响应 的 0 ~ 2 阶 谱 矩 
的 简明 闭 式 解 ,为 结构 动力 可 靠 度 分 析 提 供 了 一 种 
新 的 路 径 。 

3) 通 过 与 虚拟 激励 法 对 比 ,本 文 方法 在 分 析 结 
构 响 应 方差 和 谱 矩 时 为 封闭 解 ,具有 效率 高 和 精度 
好 的 优点 ,同时 可 用 来 验证 虚拟 激励 法 谱 矩 分 析 时 
的 精度 问题 。 

4) 由 于 本 文 复 模 态 法 对 体系 进行 解 耦 , 故 本 文 
方法 同样 适用 于 多 自由 度 非 黏 滞 阻尼 结构 基于 
Kanai-Tajimi 谱 的 地 震动 响应 分 析 。 
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